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SUSY Ward identities for the N=l SU(2) SUSY Yang-Mills theory are studied on the lattice in a non- 
perturbative numerical approach. As a result a determination of the subtracted gluino mass is obtained. 



1. Introduction 

The formulation of SUSY gauge theories on 
the lattice is problematic since the discretization 
breaks the Poincare invariance, a sector of the su- 
peralgebra. In the Wilson approach the suppres- 
sion of unphysical states in the fermionic sector 
is obtained by the introduction of an extra-term 
(Wilson term) which explicitly breaks SUSY. The 
restoration of SUSY in the continuum limit can 
be verified by considering the related lattice Ward 
identities (SUSY Wis) §. 

We focus on the N=l SU(2) SUSY Yang- 
Mills theory (SYM) (see also ||^ and references 
therein). This is the SUSY version of quantum 
gluodynamics where gluons are accompanied by 
fermionic partners (gluinos) in the same (adjoint) 
representation of the color group. As a conse- 
quence of the explicit breaking of the symmetry, 
the SUSY Wis assume in the lattice theory a pe- 
culiar form. We restrict the analysis to the on- 
shell regime . A subtracted gluino mass ms ap- 
pears; in addition, the SUSY current S^{x) gets a 
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multiplicative factor Zs and a new mixing term 
ZtO^Ti^Ix) is added to the nominal Wis of the 
continuum. 

In this contribution we present the non- 
perturbative determination of the quantities 
msZg^ and ZtZ^^ from the numerical analy- 
sis of the SUSY Wis. Preliminary results were 
presented in More details, including related 
theoretical issues, will be presented in a forthcom- 
ing publication. This study is also complemented 
by a perturbative computation |^ . 

The numerical computations were performed 
on the CRAY-T3E computers at John von Neu- 
mann Institute for Gomputing (NIG), Jiilich. We 
thank NIG and the staff at ZAM for their kind 
support. 

2. Method 

We consider the zero momentum lattice SUSY 
WI with insertion 0{y) 

^5E((^o^o(^)) 0{y))+ZTYi{V,U{x)) 0{y)) 

X X 

^ msYix\x)0{y)) + 0{a) . (1) 

This WI is valid in the on-shell regime where 
X ^ y and for gauge-invariant operators 0{x) 
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Figure 1. msZg^ as a function of t with insertion 
-^(sp)^2;) (point-split currents). 

(see § for the general case). We explain 
briefly the meaning of the various quantities 
in (0). S^ix), Tl[x) and x\x) (sink oper- 
ators) are lattice forms of the SUSY current 
S^l{x) = -(Tpff7^Tr(i^pcr(x)A(a;)), the mixing cur- 
rent T^{x) = 2'y^'TT:{F^^{x)\{x)) and the soft 
breaking operator x(x) = o'paT^T:{Fpa(x)\{x)), re- 
spectively. The trace is taken on color indices and 
A(a;) is the adjoint Majorana field of the gluino. 
We consider Q a local and a point-split defini- 
tion of the currents. The field tensor Fp^{x) is 
replaced by a clover-symmetric lattice field ten- 
sor. The quantities Zs and Zt are renormaliza- 
tions coming from the lattice SUSY breaking, ms 
is the gluino mass shifted by additive renormal- 
ization. The condition ms = is supposed to 
correspond, in the continuum limit, to the phys- 
ical situation where the gluino is massless and 
SUSY is restored. 

In our analysis we consider the lowest dimen- 
sional insertion operators Oix) {d — 7/2). One 
has essentially two choices jj]. One is the time- 
slice operator obtained from x'(^) by discarding 
time-like plaquettes, O^'^^'lx) = x'"'^^\^)] another 
possibility is 0^'^\x) = T^°'^\x), extended in the 
time-direction. We smear the insertion operators 
by combined APE and Jacobi smearing on the 
gluon and gluino fields respectively. Smearing 
significantly improves the signal for (x) but 
not for Tq°^\x). This is presumably because the 



latter contains temporal links, for which a multi- 
hit procedure is more appropriate than smearing. 
Such a procedure is however computationally too 
expensive in our setup with dynamical fermions. 

For a given insertion 0{x) the WI (Q) results 
in two independent equations when composing 
the spins of sink and insertion operators. The 
solution of the 2x2 linear system allows the 
non perturbative determination of msZg^ and 
ZxZg^ for each time-separation t = xq — yo- 
See Fig. |^ for an example. Alternatively we 
solve the overdetermined linear system for sev- 
eral time-separations (tmin, ■ ■ ■ ,Lt/2); the values 
of msZg^ and ZtZ^^ are obtained in this way 
through a least-square fit. 

3. Results 

Configurations were generated on a 12'^ x 24 
lattice at f3 — 2.3 by means of the two-step multi- 
bosonic algorithm (TSMB). See [|| and references 
therein for more details on the algorithm. See 
also Q for an application to QCD with three 
dynamical quark fiavors. The configurations at 
K = 0.1925 were produced in Results concern- 
ing K = 0.1925 and 0.194 were already presented 
in 01 . We add here more statistics at k = 0.194 
and a new simulation point, k = 0.1955. The al- 
gorithmic setup was optimized in order to reduce 
autocorrelations for light gluinos. 

In Table |l| we report the complete results for 
the global fit over a range of time-separations. 
The smallest time-separation included in the fit 
trnin was choscn such that contact terms in the 
correlations are absent; this means tmin — 3 for 
insertion x.'''^^\x) and tmin = 4 for Tq°''\x). Dis- 
cretization effects can be evaluated by comparing 
determinations from different insertions, see data 
for K = 0.1925 and k = 0.194. For k = 0.1925 
we also report results for different definitions of 
X^'*P^(a;), namely for the simple-plaquette defini- 
tion of the field tensor and for different smear- 
ing parameters. The deviation ranges between 
20% and 40% for rngZ^^. Data from T^^^^^x) 
are however subject to large statistical fluctua- 
tions and thus 0{a) effects cannot be reliably es- 
timated. 

In Fig. Ijwe report the determination of msZg^ 
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Table 1 

Summary of results. 



point split currents local currents 





operator 








ZrZg^ 


0.1925 
0.1925 
0.1925 


^(sp) (*) 

^(sp) J-^sC-J 


0.176(5) 
0.182(6) 
0.1969(47) 


-0.015(19) 
-0.044(16) 
-0.058(14) 


0.166(6) 
0.173(6) 
0.1821(47) 


0.183(14) 
0.176(14) 
0.146(11) 


0.1925 


rp{loc) 

^0 


0.132(16) 


0.11(7) 


0.144(18) 


0.29(6) 


0.194 




0.148(6) 


-0.038(19) 


0.124(6) 


0.202(15) 


0.194 


rp{loc) 

^0 


0.095(27) 


0.11(13) 


0.076(30) 


0.27(9) 


0.1955 




0.0839(4) 


-0.051(13) 


0.0532(40) 


0.179(10) 



* With plaquette field tensor. 

** With plaquette field tensor and different smearing. 
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Figure 2. msZg^ as a function of with inser- 
tion x^^^H^^) (diamonds) and Tg'"^-* (x) (triangles) 
and point-split currents. The filled circle is the re- 
sult of the extrapolation, the filled triangle is the 
determination of Kc of 0]. 

as a function of the inverse hopping parameter. 
The expectation is that mgZg^ vanishes linearly 
when K — > Kc. We see a clear decrease when k 
is increased towards Kc- An extrapolation us- 
ing data from insertion x^'^^\x) gives as a re- 
sult: Kc = 0.19750(38) for the point-split cur- 
rents and Kc — 0.19647(27) for the local ones. 
The result can be compared with the estimate 
Kc = 0.1955(5) from the study of the first order 
phase transition H. An analogous analysis for 



the quantity ZtZZ (fitting to a constant, in this 
case) gives ZtZ^ — —0.039(7) for the point-split 
currents and ZtZ^"^ = 0.185(7) for the local ones. 

Our results demonstrate the feasibility of im- 
plementing lattice SUSY Wis in order to verify 
supersymmetry restoration in a non-perturbative 
framework. 
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